Large-Scale Structure with Gravitational Waves I: Galaxy Clustering 
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Observed angular positions and redshifts of large-scale structure tracers such as galaxies are 
afTected by gravitational waves through volume distortion and magnification effects. Thus, a grav- 
itational wave background can in principle be probed through clustering statistics of large-scale 
structure. We calculate the observed angular clustering of galaxies in the presence of a gravitational 
wave background at linear order including all relativistic effects. For a scale-invariant spectrum of 
gravitational waves, the effects are most significant at the smallest multipoles (2 < < 5), but 
typically suppressed by six or more orders of magnitude with respect to scalar contributions for 
currently allowed amplitudes of the inflationary gravitational wave background. We also discuss 
the most relevant second-order terms, corresponding to the distortion of tracer correlation functions 
by gravitational waves. These provide a natural application of the approach recently developed in 
Schmidt and Jeong [J]. 



I. INTRODUCTION 

The origin of the initial perturbations which gave rise 
to the structure in the Universe is one of the most pro- 
found questions in cosmology. Currently, the most popu- 
lar scenario is inflation 0,1^, a phase of accelerating ex- 
pansion in the very early Universe which produced seed 
perturbations as quantum fluctuations frozen after exit- 
ing the horizon. One of the key predictions of inflation is 
a potentially observable background of stochastic gravi- 
tational waves (GW). A detection of a GW background 
would allow for a determination of the energy scale of 
inflation, and pose a significant challenge to competing 
scenarios for the origin of the initial perturbations. 

The polarization of the cosmic microwave background 
(GMB) is widely considered to be the most promising 
probe of the primordial GW background in the near fu- 
ture d, m. However, given the profound impact of a 
detection, it is worth studying complementary observa- 
tional techniques, in order to be able to provide inde- 
pendent confirmation of a positive result. The study of 
the effect of GW on large-scale structure observables has 
a long history. Linder Q considered the distortion of 
galaxy correlation functions and derived an upper limit 
on the GW background. Bar-Kana Q studied the ap- 
parent proper motion of distant objects induce d by GW 
(see also d, [^). More recently. Book et al. [ij have 
studied the prospects for using gravitational lensing of 
the GMB for this purpose. At redshifts of order 10 - 200, 
the 21cm HI emission from the dark ages has been pro- 
posed as a potentially extremely sensitive probe of a GW 
background [111, IS- Due to its three-dimensional na- 
ture and observable structure on much smaller scales, the 
21cm emission should in principle be able to probe GW 
amplitudes orders of magnitude smaller than the GMB. 
The shear, measured through correlations of galaxy el- 
lipticities, has also been studied as avenue for detecting 
a GW background [H, [l3|, though these authors have 
concluded that this measurement will likely not be com- 



petitive with the GMB (see also jT5|). 

The goal of this paper, and its companion [l5|, is to 
systematically and rigorously derive the GW effects on 
large-scale structure observables. While we restrict our- 
selves to a linear treatment in the tensor perturbations, 
we strive to keep the results as general as possible other- 
wise. This paper deals with observed densities of large- 
scale structure tracers, which have so far not been inves- 
tigated in the context of two-point statistics as a probe 
of GW. The companion paper deals with shear (as mea- 
sured from, e.g., galaxy ellipticities) . 

Since there is no 3-scalar that can be constructed from 
tensor perturbations at linear order without making ref- 
erence to some external (3-)vector or tensor, the "intrin- 
sic" density of tracers, i.e. the density that would be 
measured by a local comoving observer, is not affected 
by tensor modes at linear order. Thus, the impact of 
tensor perturbations is exclusively due to projection ef- 
fects which can be derived in analogy to jl6l - [l9j using 
the geodesic equation. The main observable we consider 
is the angular power spectrum C{1) of tracers. Since the 
GW effects arc most important on the very largest scales, 
the 3D power spectrum P{k) is not a meaningful quantity 
for this purpose. 

Note that in contrast, the shear is itself a tensorial 
quantity, and thus there is a possible "intrinsic" contri- 
bution correlated with the GW background, analogous 
to the intrinsic alignment effect present for scalar per- 
turbations. This issue, which has not been investigated 
before, is the topic of the companion paper [isj . 

As emphasized by Kaiser and Jaffe [2^ , there are some 
key differences in how tensor perturbations affect pho- 
ton geodesies as opposed to scalar perturbations: scalar 
modes which are transverse to the line of sight lead to 
a significantly amplified coherent deflection, whereas the 
same does not happen for tensor modes as they them- 
selves propagate at the speed of light. Furthermore, while 
scalar modes grow (in the matter-dominated era), tensor 
modes redshift away. Thus, the intuition that transverse 
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modes at lower redshifts contribute most of the projec- 
tion (lensing) effects does not hold anymore for tensor 
modes. Rather, the GW contribution to the clustering 
of LSS tracers is dominated by contributions close to the 
time of emission, and time-derivatives of the perturba- 
tions are (at least) as relevant as spatial derivatives. 

The outline of the paper is as follows. We introduce 
our notation and conventions in § |TI1 § IIIII presents the 
geodesic equation for a tensor mode and the derivation 
of the tensor mode contributions to the observed galaxy 
density, including the magnification bias contribution. 
The galaxy angular power spectrum is discussed in § IIVI 
We also highlight the differences to the scalar case in 
this calculation. § |V] deals with the relevant higher or- 
der terms neglected in § IIIII We conclude in § IVII The 
appendix contains details on some aspects of the calcu- 
lation. 



II. PRELIMINARIES 

We begin by introducing our convention for metric and 
tensor perturbations and some notation. For simplicity, 
we restrict ourselves to a spatially flat FRW background, 
and consider only tensor (gravitational wave) modes in 
the main part of the paper. The perturbed metric is then 
given by 

ds^ = a^{ij) [-dr]"^ + (% + h,j)dx'dx^] , (1) 

where hij is a metric perturbation which is transverse 
and traceless: 



(2) 



In order to simplify the analysis, we shall also consider 
the conformal metric. 



ds^ = —drj^ + {Sij + hij) dx^dx\ 



(3) 



where rj denotes the conformal time, for our analysis of 
the light deflection. 

We then decompose hij into Fourier modes of two po- 
larization states, 

/i,,(k,,y) =e+(k)/i+(k,7y) + e,f^.(k)/^"(k,77), (4) 

where efj(k), s ~ +, x, are transverse (with respect to 
k) and traceless polarization tensors normalized through 
efjC'^ = 25''^ . We assume both polarizations to be 
independent and to have equal power spectra: 

{hs{k, r^)h,, (k', r/)) = (2^)35^(k - k')6ss' ^Prik, ry, r?')- 

(5) 

Here, rj denotes conformal time, and the unequal-time 
power spectrum is given by 

PTik,v,v')=TTik,i^)TTik,i^')PToik), (6) 



where TT{k,r]) is the tensor transfer function, and the 
primordial tensor power spectrum is specified through 
an amplitude and an index ht via 



PTo(fc)=27r^fc-^( A 



(7) 



Following WMAP convention [21[, we choose kg ~ 
0.002 Mpc^^ as pivot scale. Throughout, we will assume 
a scalar-to-tensor ratio of r = 0.2 at fco (consistent with 
the 95% confidence level WMAP bound), which together 
with our fiducial cosmology determines A^. The tensor 
index is chosen to follow the inflationary consistency re- 
lation, Tlx = —r/8 = —0.0025. For the expansion history, 
we assume a flat ACDM cosmology with h = 0.72 and 
ilm — 0.28. Contributions from scalar perturbations are 
evaluated using a spectral index of = 0.958 and power 
spectrum normalization at z = of cts = 0.8. 
From Eq. (|4|) and Eq. ([5]), we easily obtain 



(/i,,(k, 7j)hki{k', Tj')) = (27r)3<5^(k - k') 



(8) 



e+(k)e+(k) + e><.(k)e,^,(k) 



X ■^PT{k,T],v') 

{h,,{k, v)h'Hk', 77')) = {2TTfSD{k - k')PT(fc, V, V')- 

Long after recombination, the transverse anisotropic 
stress which sources gravitational waves becomes neg- 
ligible, and the tensor modes propagate as free waves. 
During matter-domination, the tensor transfer function 
then simply becomes 



TT{k,v) 



(9) 



which however is still valid to a high degree of accuracy 
during the current epoch of acceleration. We will use 
Eq. ^ throughout. 



III. TENSOR CONTRIBUTIONS TO THE 
OBSERVED GALAXY DENSITY 

In this section we derive the GW contribution to the 
observed density of tracers, including the magnification 
bias effect. We follow the notation in Jeong et al. [l^ 
(see also [l|) which is summarized in their Sec. II A. The 
zcro-th order photon geodesic in conformal coordinates 
[Eq. (jS])] is simply 



2;^(x) = iVo - X,nx) , 



(10) 



where the comoving distance x along the geodesic serves 
as affine parameter, with x = corresponding to the 
observer's location. Here and throughout, n denotes the 
unit vector in the direction of the observed position of 
the source {n in [3). Hence, 



dx 



(-l,n). 



(11) 
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In the following, z will stand for the observed redshift of 
the source, and x = x(^) is the conformal distance corre- 
sponding to that redshift when evaluating the distance- 
redshift relation "xiz) in the background. 

We decompose vectors into transverse and longitudinal 
parts with respect to the line-of-sight. 



X 1 = X'' - 



(12) 



and correspondingly define longitudinal and transverse 
derivatives through 



(9|| = fi^ di 



(13) 



(14) 



(15) 



Finally, we define 
and make use of 

diifi' =fl'^±^ = 0. 
A. Photon geodesies with a tensor mode 



Wc now briefly outline the derivation of the displace- 
ments Ax* of the true emission point from the obscrva- 
tionally inferred position. This is a special case of the 
derivation in App. B of Schmidt and Jeong [l[, to which 
the reader is referred for more details. We parameterize 
the linear order deviation of the photon geodesic as 



—— = -1 + (5i/,n 
dx 



(5e) 



(16) 



The initial conditions for integrating the geodesic equa- 
tion are set by demanding that the components of the 
photon momentum measured in a locally orthonormal 
frame at the observer's location match n. For that, we 
construct an orthonormal tetrad (e")^ (a = 0, 1, 2, 3) car- 
ried by an observer so that at the observer's location 



(17) 



Then, the photon four-momentum measured by the ob- 
server is given by 



(l,nO-((e°)^,(e%) 



dx ' 



(18) 



A detailed calculation is presented in appendix B of [l|. 
Using that the four-velocity of comoving observers is 
given by m'' = a~^(l, 0, 0, 0), this leads to 



Su{x^O)^0 
6e\x = 0)^ -l{h^,)on\ 



(19) 



Here and throughout, a subscript o indicates that the 
quantity is evaluated at the observer's location. The cor- 
responding initial condition for the geodesic equation for 
comoving observers with general metric is presented in 
Eq. (Bll) of Q. In case of scalar perturbations, includ- 
ing (5e'(x = 0) is important to ensure gauge- invariant 
expressions [l9|. While there is no gauge ambiguity in 
tensor modes, we will show that the observer term is nu- 
merically important for the quadrupole of the observed 
galaxy density. The redshift perturbation is related to 
the shift in the frequency through 



Sz{x) = —Siy = 

where we have defined 

/ill = hi 



1 



Kdx, 



(20) 



(21) 



Here and hereafter, a prime denotes a derivative with 
respect to conformal time, if not used for distinguish- 
ing different variables. The distinction between the two 
should be clear from the context. The redshift of the pho- 
ton along the perturbed geodesic is given by 1 + z{x) = 
[l + 6z{x)]/a{x^{x))- Requiring that the redshift at emis- 
sion equals z yields an equation for the first-order per- 
turbation to the affine parameter at emission Xe p^ . 



Xe=X + Sx 

Sx = Sx° - \^Sz, 
H{z) 



(22) 



where 5x^ is the perturbation to the 0th component of 
the geodesic evaluated at x = X- Wc now relate the 
observed position x, inferred assuming the unperturbed 
geodesic x^, and the true position x through (see Fig. 1 
in [3) 



Ax = X — X = 5x[x) + n(5x, 



(23) 



where 5x{x) is the spatial perturbation to the geodesic 
evaluated at x = X- We can then decompose the dis- 
placement Ax into perpendicular and longitudinal parts, 

A.T|| = 5x'n, + 5x° - ^^7^5z (24) 
ti yz) 

Aa;^ = Ax* -n^Axii. (25) 



Explicitly, 



and 



dxh\\ 



l + z 
2H{z) 



dxh\ (26) 



Axl=-x[(/i.,)>^"-(/i||)„^ 



dx 



ai/i|i + ^(/i||n'-/iyn^) 
^ X 



(27) 



4 



Note that Eqs. (|26)) - ([27)) can be obtained directly from 
Eqs. (43)-(47) in [l|~~(with 5z = A In a) by setting 
A = _B = w = 0, or, alternatively, from the correspond- 
ing Eqs. (38)-(39) in [l^ by restricting to a transverse- 
traceless metric perturbation and setting = ft.||/2. 

In the following, we will further need the convergence 
k, defined through (see App. 



1 



dx 



3, 

X 



-V2 



ax— ^11- 

XX 



(28) 



Here Vj^ = X^'^\ is the Laplacian on the sphere. Also, 
we will use 



1, l + z 

-"II 

2 " 2H{z) 



h'n 



H{z) d 

2 ai 



l + z 



H{z) 



dxK, (29) 



where we have used dx = dz/H{z). While k is the usual 
coordinate convergence, 9;^Aa:;|| is the distortion of the 
volume along the line of sight and can thus be seen as a 
"radial convergence" . 



B. Observed galaxy density 



The observed comoving number density of galaxies 
a^hg is related to the true comoving number density a^rig 



through 



a {z)ng{x, z) 



a {z)ng{x, z) 1 + 



dAx' 



(30) 

The first term in brackets comes from the covariant 
volume factor \/\g\ and is equal to unity, since for 
transverse-traceless metric perturbations Sg^^^ = at lin- 
ear order. The factor a^{z)ng{x,z) is the true comoving 
number density at the point of emission, which we ex- 
pand as 



a^(z)ng(x, z) = a^{z)ng{z) [I + (5g(x, z)], 



(31) 



by defining the intrinsic perturbations to the comoving 
number density Sg. z is the redshift that would be mea- 
sured for the source in an unperturbed universe, and is 
related to z through 



l + S = {l + z){l + Sz). 



(32) 



Note that, when inserting Eq. ([31]) into Eq. ([5(7)) . the dis- 
tinction between x and x in the argument of Sg is second 
order if we regard intrinsic galaxy density perturbations 
as first order, which we will do in this section. The rel- 
evant additional terms will be studied in § [Vj Finally, 



1 -I- 9iAx* is the volume distortion due to gravitational 
waves, which as derived in Il9| becomes 



dx^ 



d^Ax 



2Ax 



2k. 



X 



(33) 



Thus, gravitational waves affect the observed density of 
galaxies through a volume distortion effect, and by per- 
turbing their redshifts so that we compare the measured 
galaxy density hg to the "wrong" background density 
ng{z). The latter effect is quantified by the parameter 



_ dln(a3 



din a 



-(1 + 5) 



,.dln(a3 



dz 



(34) 



Note that this parameter can be measured for a given 
galaxy sample, provided the redshift-dependence of the 
selection function is understood. 

We can now summarize the tensor contributions to the 
observed galaxy density perturbation as 



6gT{x, z) = beSz — 2k — z 
X 



^ 1 + z 

'^''^^-'imJo '^''i 



2H{2) II 2 82 



l + z 



Hiz) 



dxh\y 

(35) 



Here, the subscript T denotes tensor contributions, and 
we have assumed that the intrinsic density perturbation 
5g does not correlate with the tensor modes (following 
the arguments in § |l| . Note that hij only enters through 
This has to be the case, since the galaxy density 
is a scalar quantity and h\\ is the only non-trivial scalar 
linear in hij . The latter also implies that tensor modes do 
not contribute to the monopole and dipole of the galaxy 
density. 



C. Magnification bias 

In the last section, we have assumed that all galax- 
ies are included in the sample. In reality, most large- 
volume surveys are limited in flux. A cut on observed 
flux induces additional fluctuations in the galaxy density, 
since perturbations to the photon geodesic (e.g., gravita- 
tional lensing) modify the observed flux of a given source. 
The magnification A4. can be derived as the perturba- 
tion to the angular diameter distance squared ( [l9[ , with 

<5A4 there = Ashore), 



l + M 



D2 



(36) 



where Da is the true angular diameter distance to the 
source while Da{z) is the background angular diameter 
distance-redshift relation. Alternatively, one can use the 
standard ruler approach of [l[ , which for a metric of the 
form Eq. ([T|) yields (see Eq. (105) in [l[, where we set 
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d In To /d In a = neglecting the higher order effect of an 
evolving number count slope) 

Xt = -2(5z + -/i|| ^+2k. (37) 

2 X 

We then parametrize the effect on the observed galaxy 
density through a parameter Q, 



^gT ~ '5gT(no magn.) + QMt- 



(38) 



For a purely flux-limited survey, Q = — dlnng/dln/cut, 
where /cut is the flux cut. More generally, Q can also 
receive a contribution from a size cut 1231. 



D. Summary of tensor contributions 

Combining the result of the last two sections, we ob- 
tain the following expression for the linear-order tensor 
contributions to the observed galaxy density: 



IgT = (be - 2Q)6z - 2(1 - Q)>i - ^— 



l-Q 
X 

H{S) d 

2 9i 



S 1 + z 

'^^'ll + ^io 



1 + S 



H{~z) 



dxh'ii ■ 



2H{z) II 



(39) 



For later convenience, we reorder the terms as follows: 

4t = hh + m + f J '^h + fj jh (40) 

Here, all terms outside integrals (without subscript o) 
are evaluated at %) and the integrals go from to x- The 
coefficients are given by 



f'~ 

•I V 



i(Q-i) 

l+~z 



(41) 



2H 

f = Q-i 

1 



l-2Q+(l + z)^)-(Q-l) 



H 



l + z 



Hx 



^(2-1) 



/o = 2(2-1)- 



IV. OBSERVED GALAXY POWER SPECTRUM 

A. Angular power spectrum 

Consider a galaxy sample with a redshift distribution 
dN/dz, normalized to unity in redshift. Then, the pro- 
jected galaxy overdensity as a function of position on the 
sky is given by 



r 

Ag(n) = y dz — Sg{x{z)h;z), 



.dN~ 
dz 



(42) 



We will assume that the quantities b^, Q describing the 
galaxy sample are independent of redshift for simplicity. 
We can then write the multipole coefficients of the galaxy 
density as 

«L= J d^nYr„,{h)Agin). (43) 
We can write all individual contributions to Eq. ((5^ as 



^(n,x)=: / dxWAix,x)h\\ixi^,x) 
Jo 



dxWA{x,x) 



d^k 
(2^ 



e'''-^^h'h^h,j{k,fio~x)- (44) 



Note that terms involving can be brought into the 
form A{h) by including d\nTT{k^'q) / drj in Wa{x)- The 
observer term 5/i||o/3 contained in k can similarly be 
written with Wa{x) ~ 5/3(5_d(x). We will deal with 
that term specifically in § IIVBI By changing the order 
of integration, we can then write the contribution to the 
projected galaxy overdensity as 



r°° dN 

Wa{x)= / d~z—WA{x,xm^ 
Jzix) dz 



rfx WA(x)^||(xn,x) 



(45) 



Note that ff Wa = Sd{x - x), Wa(x) = iHdN/dz)Ui^^y 
We now consider the contribution of a single plane wave 
tensor perturbation with k-vector aligned with the z- 
direction. Then, 

n^n^hij{'k,r]) = sin^ 6* [cos 20 /i+(k, 7?) -f sin 2(?!) /i^ (k, 7?)] 



sm t/ e 



where 



/il,2 =-(/!+ ±l/l^). 



(46) 



(47) 



Note that the power spectra of these circular polarization 
staters are P,,^,,^ = Ph^h^ = -Pt/8, while Phih2 = 0. Let 
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us denote as ^(n,k) the contribution to A{ii) from this 
plane-wave tensor perturbation. We have 

A(n,k) = J dxWA{x)e"'''^il-^^^) (48) 

X [e^'^h,{k, Vo-x) + e-2'*/i2(k, Vo - x)] 

where n = cos 9 is the cosine of the angle between fi and 
k. Note the e*^"^ factors which are the key difference 



to the case of scalar perturbations. The multipolc coeffi- 
cients of A are then obtained as follows: 



afm^ J -^af^i^), (49) 



at{k) = I d2fil^:„(n)A(n) 

= J dxWAix) J <fhY:^it^,^)e'''^'^{l~fi') P*/ii(k,7yo-x) + e-''^/i2(k,r;o-x)] . (50) 
We now use the relation (see App. Al in [l[) 

I dnYr^{l-^^')e^^'^e^^^=- ^4^(21 + 1),|^±^^^^J^S^±,, (51) 

which yields 

aL(k) = -J^^l±^{4n)^' J dx^Aix) [hi{Km - x)S,n2 + h^iKv^ ~ x)Sm-2]-^^. (52) 



Further, we can use the properties of the spherical har- 
monics to obtain 



ai, 



(53) 



which we will use to evaluate the convergence contribu- 
tion. Using the definition of the angular power spectrum, 
we can now easily write down the cross-correlation be- 
tween two different projections A, B of h\\ (App. Al in 
i): 



AB _ 



21 
1 (/ 



2)! 



(k) 



271 {I - 2)! 

dx^x{x)TT{k,r]o 



k^dkPTo{k)F/^{k)Fi^{k) (54) 



X) 



jijkx) 
{kxY ■ 



(55) 



It is then straightforward to evaluate auto- and cross- 
correlations of the observed angular correlation of galax- 
ies, neglecting the contribution of the term for the 
moment (see § IIVBp . In particular, if we want to evaluate 
the total tensor contribution to the angular (auto-)power 
spectrum of galaxies, we set Wa = Wb 

d In Tt 



Wg, where 



Wg{x)=hSD{x-x) + 



drj 



f- + /4 - U lii + 1) 

X X 



no-x 
X - 



{f~SD{X'X) + f') 



XX 



(56) 



and the coefficients are defined in Eq. pT|) . Note that 
the divergent pieces oc x~^ cancel for 1^2: 



f - 6/k = 0, 



(57) 



while for Z > 3 the Bessel function in Eq. ([55]) ensures 
that all terms are regular. Fig. [1] shows numerical results 
for a galaxy sample with a sharp (observed) redshift of 
z ~ 2, and for be ~ 2.5, Q = 1.5. The colored lines 
indicate the separate contributions proportional to pro- 
jections of and V^/i||. While for Z < 4 all terms 
contribute significantly, the contribution dominates 
at higher I. This contribution is the same as the GW 
effect explored in Q . We will return to this in § IIVDI 
Note that the total contribution (black solid line) is sig- 
nificantly smaller than the individual contributions for 
I = 2, 3. We will discuss this in the next section. 



B. Quadrupole 

Fig. [T] shows that the individual tensor contributions 
to the observed angular power spectrum increase rapidly 
towards small I, while the total contribution is much 
smaller. These cancellations are essentially a conse- 
quence of causality, which demands that the observed 
clustering of galaxies cannot depend on tensor modes 
that are super-horizon today. 
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FIG. 1: Contributions to the observed galaxy angular power 
spectrum from inflationary gravitational waves, for a sharp 
source galaxy redshift of 5 = 2, and using the tensor mode 
power spectrum defined in § [Til The black solid line shows the 
total contribution, while the colored lines show contributions 
proportional to line-of-sight integrals of /ly (blue dotted), fty 
(green short-dashed) and Vq/i|| (magenta long-dashed). Here, 
we have assumed be = 2.5, Q = 1.5. The black square at I — 2 
indicates the result for Z = 2 if the observer term is neglected 
(see ^ HVB)) . 



In particular, it is important to take into account the 
last term in Eq. which quantifies the shearing of the 
observer's coordinate system by the gravitational waves. 
We refer to this as the observer term. Since hij{o) is 
a constant (transvcrse-traceless) tensor, this term only 
contributes to the quadrupole I ~ 2. Specifically, we 
have an additional contribution to Eq. (|55|) for galaxies. 



Fr(fc)=/olimrT(fc,%-x) 



{kx? 



= Y^foTT{k,i-io)- 



(58) 



The solid black line in Fig. [T] includes this term, while 
the black square indicates the value of Ci=2 we would 
obtain without this term. Neglecting the observer term 
results in an overestimation of the tensor contribution 
to the galaxy quadrupole by three orders of magnitude. 
The reason for this significant effect becomes clear when 
considering the contributions to Fj^^^ ^ function of k 
(Fig. [2]). The individual contributions to Ff' approach 
a constant as fc ^> 0, while the sum goes to zero for 
k/Hg < 1, as demanded by causality. When neglect- 



FIG. 2: Contributions to the kernel Fi(k) for Z = 2 (thick) and 
/ — 20 (thin, only total contribution shown, scaled by 10*), 
for a sharp source redshift 5 = 2 and the same parameters 
as in Fig. [T] Note that the separate contributions have non- 
zero weight for fc — > 0, while the total Fj^ (black solid) is only 
non-zero for k > 10^* ft/Mpc, as required by causality. 



ing the observer term (light blue in Fig. [2]) on the other 
hand, a residual constant contribution to i^;^2 remains 
for fc — > 0, which together with the steeply falling tensor 
power spectrum leads to a significant overestimation of 
the quadrupole. 



C. Limber approximation 

In the context of angular galaxy clustering, one often 
uses the Limber approximation [23| which significantly 
simplifies the calculation of Cf . The underlying assump- 
tion is that the dominant contribution to the angular 
clustering comes from galaxy pairs that are at similar dis- 
tances along the line of sight. It is instructive to consider 
this approximation in the context of tensor modes. Since 
the Limber approximation works best for a broad redshift 
distribution, we will here consider a redshift distribution 
roughly as ex pec ted for the Large Synoptic Survey Tele- 
scope (LSST [H), 



dN 
dz 



oc z exp 



(59) 



with zq ~ 0.15 and /3 = 0.73, yielding a mean redshift of 
1.2. 

The Limber approximation can formally be applied by 
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using 



edkF{k,x)ji{kx)ji{kx') 



1+1/2 



X 



Sd{x-x')- 
(60) 

In the usual application to scalars, the funetions of fc, x 
involved (apart from the Bessel functions) are smooth 
and positive, whereas for tensor modes, the transfer func- 
tion is oscillatory. Performing the k and then one of the 
X integrals in Eq. ([54|) leads to 



AB 



{1 + 2)1 



1 



1 



{1-2)1 {I + 1/2)^ A 
xWa{x)Wb{x) 



dx 

,,2 



Pto[^-^] (61) 



Tt {1 + 1/2) 



Vq-X 
X 



1 2 



Here we have used that Tt is only a function of kr/. Note 
that for ^ :§> 1, the prefactor approaches 1/4. Given 
that Pro oc we immediately see that the Limber 

approximation predicts Cf oc ;~3+nT j-qj. \g^Ygc I, if Wa.b 
are /-independent. 

Fig. ini shows the different contributions to Cf (as in 
Fig. [T]) for the rcdshift distribution Eq. (|59p using the 
full calculation (thick lines) and using the Limber approx- 
imation (thin lines). Clearly, the Limber approximation 
predicts the wrong Z-scaling of all terms, and we do not 
see any improvement in the approximation for high I as 
in the scalar case. Thus, the Limber approximation is not 
applicable for calculating the tensor contribution to the 
angular galaxy power spectrum at any £. There are two 
reasons for this. First, the tensor mode power spectrum 
is falling as ^ fc~'^, so that the assumption that pairs of 
comparable line-of-sight distance dominate because they 
have the smallest separation does not hold. Second, ten- 
sor modes oscillate and decay towards late times (while 
scalar modes grow), so that the contributions to Cf are 
concentrated at large scales and high rcdshifts near the 
source. 



D. Dependence on galaxy sample 

Fig. d] shows the tensor contribution to Cf for differ- 
ent source redshifts. Here, we have assumed a Gaussian 
redshift distribution centered on z with RMS width of 
0.03(1 -|- z), emulating the effect of photometric redshift 
errors. We have kept be ~ 2.5, Q = 1.5 fixed. While 
the low multipoles {I < 4) do not depend very strongly 
on z, the contribution for higher / grows by an order of 
magnitude when going from z = 1 to z = 2, and from 
z = 2 to z = 5. We found that reducing the width Az 
of the redshift distribution has no impact on the results 
at relevant I. The reason is that the bulk of the tensor 
mode contribution comes from near horizon-scale modes, 
as tensor modes decay once they enter the horizon. Such 
long- wavelength modes arc not affected by averaging over 
a reasonably narrow rcdshift window (Az < 1). 
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FIG. 3: Contributions of tensor modes to the angular galaxy 
power spectrum for a broad redshift distribution (Eq. (|59p ) 
expected for LSST, separated into different contributions as in 
Fig. [T] (again using be = 2.5, Q = 1.5). The thick lines show 
the exact calculation, while the thin lines show the Limber 
approximation using Eq. (|61fl . 



Fig. [5] shows the effect of varying the galaxy sample 
parameters b^, quantifying the redshift evolution of the 
average number density, and Q, which determines the 
magnification bias contribution. We see a significant ef- 
fect at low multipoles when varying Q and especially be- 
This is not surprising since we have seen that there is a 
high degree of cancelation between different terms at low 
I (Figs. [1] and [5]), so that varying the coefficients of the 
different terms can have a large impact. 

On the other hand, is insensitive to changes in be 
and Q for / > 10. While surprising initially, this fact can 
be understood straightforwardly. At high I, the contri- 
butions from /ijj dominate (Fig.[T]). Specifically, following 
Eq. (|39p there are two contributions, f^h'^^ and /' / /iydx- 
The former contribution is a pure line-of-sight volume 
distortion effect, and the coefficient oc H~^{z) is in- 
dependent of the galaxy sample. Now consider a single 
Fourier mode contributing to both terms. Neglecting fac- 
tors of order unity, we can estimate the ratio between the 
two terms as 



rsih\dx 

m 



X))dx 



T!r{k{m - X)) 



(62) 



where 



dTT/drj. Since decays rapidly once the 
argument kr] become of order unity, we can approximate 
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the numerator as 



Ti^iHm - X))dx - T^iHvo - x))Ax, (63) 



where ~ A{kri)/k ^ 1/fc is the range in x that con- 
tributes to the integral. Thus, we arrive at an order-of- 
magnitude estimate of 



flih\^dx H{z) 



f h' 



(64) 



Physically, the integral of along the line of sight leads 
to significant cancelations that become more severe as 
k increases. It is not surprising then that the volume 
distortion term dominates at high /. Specifically, 

using the kernel Ff{k) for / = 20 shown in Fig. [21 we can 
estimate that the typical wavenumber of tensor modes 
contributing to multipole / is 



hyp ~ 0.003 h/Mpc X 



I 

10' 



(65) 



Eq. ([M)) then says that we expect to dominate for 

I > 10, in good agrecemcnt with Fig. [5] 

Thus, the line-of-sight volume distortion c>c h'^^/H is 
the single dominating term in at high Z, which ex- 
plains why the tensor contribution is independent of the 
galaxy sample at high multipoles. Note however, that the 
prospects for detecting this effect are mostly restricted to 
I < 10 anyway due to the steep decline of the signal. 



V. HIGHER ORDER TERMS 

So far. we have only kept terms linear in all pertur- 
bations. While all terms second and higher order in hij 
are likely irrelevant given the small amplitude of tensor 
perturbations, there are terms of order hij Sg which can 
be much less suppressed on small scales where Sg can be- 



come order unity. There are more terms of order hij A 
where A stands for any scalar metric perturbation. How- 
ever, those terms are much smaller than hij Sg and can be 
neglected. Note that as long as scalar and tensor modes 
do not correlate, the lowest non-trivial statistic induced 
by these terms is a four-point function oi Sg. A specific 
quadratic estimator for the tensor modes can be con- 
structed based on these terms [ill ■ Physically, one 
uses the anisotropy of the small-scale correlation func- 
tion of Sg to search for coherent large-scale distortions 
induced by tensor modes. Most of the signal-to-noise for 
detecting a GW background is contained on very small 
scales 

Noting that x = x - Ax and 1 + z = (1 + z)(l - Sz), 
the relevant terms neglected in going from Eq. (j30p to 
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FIG. 4: Total tensor mode contribution to the angular galaxy 
power spectrum as a function of source redshift, for a Gaus- 
sian redshift distribution centered on z with an RMS width 
of 0.03(1 + I). Here, be = 2.5, Q = 1.5 fixed. 



Eq. pS]) are given by 

SgT.2nd =Sg[x- Ax, (1 + z)(l - Sz) -1]-Sg (x, z) 



= — Ax • V(5g(x, 2) — Sz 



91n(l + z) ^^^"'')' 



(66) 



where Sg is the intrinsic fractional perturbation in the 
comoving galaxy density. Note that both Ax and Sz are 
linear in hij. If we approximate V(5g Sg/r, where r is 
the scale on which the correlation function of Sg is mea- 
sured, then the second term is suppressed with respect to 
the first term by a factor of r/x, since dhiSg/d\n{l + z) 
is typically order unity. This term is thus highly sup- 
pressed on small scales where the second-order terms be- 
come relevant. While the intrinsic two-point function 
^(r) of the tracer is isotropic and location- independent 
(neglecting redshift-space, tidal distortions and intrinsic 
anisotropy /inhomogeneity [TTI . [25j). the observed corre- 
lation function 



|(n,2~;n',2~') = (^(i)^(i')) 



(67) 



is anisotropic and depends on the location x. This can be 
used to measure the distortions by tensor perturbations. 
Specifically, since we have 6 free parameters in x, x', we 
can measure 6 components of the distortion field. At this 
point, it is important to stress that the terms in Eq. (|66p 
are not observable directly (as a simple example, consider 
the case of a constant deflection Ax = const). In order 
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FIG. 5: Relative impact on the total tensor mode contribution 
to galaxy clustering when changing fee from 2.5 to 4 and 1 
(thick and thin green long-dashed, respectively), and when 
changing Q from 1.5 to 3 and (thick and thin red dot- 
dashed, respectively). Here, a Gaussian redshift distribution 
centered on 5 = 2 with RMS width 0.03(1 -I- z) was assumed. 



to determine which quantities are actually observable, 
consider contours of constant ^, 

|(n,z;n,z') =^0- (68) 

These contours correpond to a fixed physical scale ro 
(on a constant-proper-time hypersurface) at the source 
through 

^(^o) = $0. (69) 

In other words, the intrinsic homogeneous and isotropic 
correlation function ^ (r) is supplying us with a "standard 
ruler" tq (or, a set of standard rulers as we are free to 
vary ^o)- In Schmidt and Jeong [l|, we carefully define 
a general standard ruler and derive the properties of the 
deflection field which are observable through it. As we 
have seen here, the distortion of correlation functions by 
tensor modes is one application of the results of [l[ . 

Finally, we point out that a non-zero three-leg coupling 
{5{ki)5{k2)hij{kz)) present at an early stage of the uni- 
verse can also imprint its signature as a local departure 
from statistical homogeneity. The optimal estimator for 
the amplitude of tensor perturbations given such a cou- 
pling has been constructed in (25| . 



FIG. 6: Comparison between tensor mode contributions to 
the galaxy power spectrum (black solid) for z = 2 (cf. Fig. 3]), 
and scalar contributions for a linear bias b — 2 (see App. |B] 
for details on the calculation of the scalar contributions). 
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FIG. 7: Comparison between tensor mode and scalar contri- 
butions to the angular cross-correlation between two widely 
separated redshift bins (5 = 1, 5' = 4) and for Q = 1 so that 
most magnification contributions drop out (linear bias = 2; 
other parameters as in Fig. [4]). 
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VI. DISCUSSION 

We have derived the complete tensor eontributions to 
the observed galaxy density at linear order. The re- 
sult is summarized in Eq. p9p . At this order, gravi- 
tational waves do not perturb the intrinsic physical den- 
sity of tracers; thus all contributions are projection ef- 
fects from the effects of GW on the propagation of light. 
We have found that, contrary to gravitational lensing 
by scalar perturbations, tensor perturbations contribute 
mainly at redshifts close to the source redshift. Together 
with the scale-invariant power spectrum of GW, this re- 
sults in a steeply falling angular power spectrum of the 
tensor contributions, with multipoles I 10 already be- 
ing suppressed by an order of magnitude with respect to 

Fig. ini shows a comparison of the tensor contributions 
with the scalar contributions to the galaxy density. Here, 
we have assumed a linear bias of 5 = 2, and all relativis- 
tic corrections are included following [l^ (see App. |B] 
for details). Clearly, the tensor contributions are sup- 
pressed by ~7 orders of magnitude with respect to the 
scalar contributions at the largest scales, for the max- 
imum currently allowed value of r = 0.2. One might 
wonder whether galaxy cross-correlations, i.e. between 
different redshift bins, could be more promising. Af- 
ter all, when cross-correlating widely separated redshift 
slices, there is little intrinsic correlation of galaxies, and 
the main contribution comes from lensing (magnification 
bias effect oc 2(Q— 1)k). Hence, the most optimistic case 
for detecting tensor modes through their modulation of 
the galaxy density would consist of cross-correlating two 
galaxy populations widely separated in redshift, both of 
which have Q = 1 so that the magnification bias effect 
drops out (this could be achieved, for example, by select- 
ing galaxies on surface brightness). This most optmistic 
case is shown in Fig. [71 We find that the tensor contri- 
bution can become as large as 10"'^ times the residual 
scalar contribution; however, this is still much too small 
to be detectable. 

Thus, given that we are not able to directly distinguish 
between scalar and tensor contributions to the galaxy 
angular power spectrum, we do not expect much detec- 
tion potential for gravitational waves from the leading 
order effect on the angular power spectrum of galax- 
ies. However, there are terms of order hij Sg (where Sg 
is the intrinsic galaxy overdensity) which induce a par- 
ticular four-point correlation function in the observed 
galaxy density Sg. This is the effect exploited in [TTI. [T2j 
for projected constraints from 21cm emission from the 
dark ages. In essence, the intrinsic (homogeneous and 
isotropic) correlation function provides us with a stan- 
dard ruler that allows us to observe certain properties 
of the distortion field, as derived in [l|. Thus, one can 
apply the the scalar- vector-tensor decomposition on the 
sky described in [l|. Most importantly, both the (2- 
) vector and (2-) tensor components allow for an E/B- 
decomposition, so that any scalar contributions to the 



distortion (e.g. from lensing or redshift-space distortions) 
do not contribute to the S-mode at linear order. We 
leave a more detailed investigation of these possibilities 
for future work. 



Acknowledgments 

We would like to thank Chris Hirata and Marc 
Kamionkowski for helpful discussions. FS is supported 
by the Gordon and Betty Moore Foundation at Caltech. 
DJ was supported by NASA NNX12AE86G. 



Appendix A: Convergence 

We define the convergence through the divergence of 
the transverse displacements: 



(Al) 



First, taking the transverse divergence of Eq. (|27)) yields 



-2k = -X {h'jd±^n^ 



hudj_^n') (A2) 



Jo L 

+ ^(x-x)|vi/.|| 

Using the fact that hij is transverse and that = 9|| 
drj , straightforward algebra then yields 



k — —h\\ 
4 " 



-/ill 



dx 



+ -^11 

A 



4^ Jo XX " 



(A3) 



The last term is familiar as the one dominating on small 
scales (Newtonian limit) for the scalar case, and is shown 
separately as blue dotted lines in Fig. [1] and Fig.[3l Note 
that due to the very different power spectrum of tensor 
modes, and the fact that tensor modes decay with de- 
creasing redshift, this term actually becomes suppressed 
on small scales. 



Appendix B: Angular galaxy power spectrum from 
scalar perturbations 

In this section, we provide expressions for the angular 
power spectrum of galaxies due to scalar perturbations 
(see Fig. [5] and [71), including all relativistic corrections as 
derived in |16l-ll9l| . 

It is convenient to use expressions in conformal- 
Newtonian gauge, where we write the metric as 

ds^ = a^{rj) [-(1 + 2^')dr]^ + (1 + 2^)6ijdx'dx^] . (Bl) 



12 



Further, for convenience we write the vefocity Vi in 
terms of a scalar velocity potential V, Vi = V^i, and re- 
late y, (f>, to the density contrast in synchronous- 
comoving gauge through (see [26l.[27|) 



4'(k,77) - *(k,7?) 



aHfk-^D{air,))5'^{k,T^o) 



4>(k, 77) + vl/(k, 77) = gik, ri)D^_ {k, 77)5^^(k, 770) 



l]i?*_C(k,?7o) 



vI/(k,77) = -[.g-l]i?<,_C(k,77o) 
($ - vl/)'(k, r,) = Asw(fc, r/)<5^^(k, 7,0), (B2) 

where / = din D/ din a, D{a) is the matter growth fac- 
tor (normalized to unity at a — 1) and we have defined 
general coefficient functions D$_, g, Disw to allow for 
non-standard cosmologies. In a ACDM cosmology, we 
have 



^*-(k,77) 



0. 



(B3) 



where 



5z=d\\V 



7^ / rfx(x-x)^vi($-vi/) 

2 Jo X 



(B4) 



denote the redshift perturbation and convergence, re- 
spectively, in conformal-Newtonian gauge. Assuming a 
linear bias relation in synchronous-comoving gauge (as 
discussed in [l9j). and expanding the different contribu- 
tions, we obtain 



Here, the subscript denotes that the quantity is eval- 
uated at z ~ 0. Throughout, we will drop observer 
terms corresponding to the monopolc and dipole, but 
keep terms that contribute to the quadrupole. 

Converting Eqs. (55) and (72) in [l^ from 
synchronous-comoving gauge to conformal-Newtonian 
gauge, using the relations given in App. D of [l^, we 
obtain (see also [13]) 



h,{5z - aHV) - -^df^V 



- 1- 
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2($ 

X 
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2$ 



1 dH 

aH dz 
-aHV) 



aHx 



[Sz ~ aHV] 
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- / dx(* - *) - 2k 
X Jo 



($ 

- 2aHV 



aHV)' 
2 
X 



+ QM 
V + 2k 



2 

+ - ($ - ^')dx + 2 
X Jo 



1 



aHx 



- 1 



{Sz - aHV), 



6e-l-2Q+47^ + 2(S-l)^ 
aH az o-Hx 



aH 



(3 - b,,)aH 



dH 
dz 



V^ + 2(l - Q)$ 



2(Q - 1) r ^($ - *) + 2(Q - 1)k 
Jo X 



df,V + i<^> + aHVy 



(B5) 



We can now evaluate the contribution of a sin- 
gle plane wave along the z-axis, i.e. (5^j(x, 77) = 
Sj^^{'k,riQ)D{a{r]))e^^^. Further, we write x = nx, and 
X ■ h ~ xfi with X = kxi x = kx- First, we have 
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aHfP ^ 1 
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dxDiswe 



'IX fl 



K = -l{l + 1) / dx^-^D^_ {k, r^o - X)e"^ C(k, ?7o) 
dlV^{aHfD),dle^'^6l^,(k,rjo) 



{fDh 



dlnaH dlnfD 
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dim 
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Following the procedure described in App. Al of [l[, it is then straightforward to derive the angular power spectrum 
of Sg, for a sharp source redshift z, in terms of the matter power spectrum today Pm(k) in synchronous-comoving 
gauge: 



k^dkP„,{k)\Ff{k)\ 



F,3ik)^bDi~a)Mx) + 



(3 - &e) 



^ ^ dH ,^ , 1 
^^-1-22+^^ + 2(2-1)-^ 
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1 dH 
aH dz 
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{fD),ji{i) + 2(1 - Q)-([,9 + Ji(i) 
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Jo XX 
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aHfD 



X • ■ 2 
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aH 



-ifD),dlji{x) 
dlnaH dlnfD' 



dim 



din a 



^ ds - 2^9 -'^)D^_ ] ji{x)+ \ dxDiswjiix), 



(B7) 
(B8) 



where again x = kx, x = kx. We can also easily derive the power spectrum of the magnification itself, which is an 
observable (see [l|): 



Cm{1) 



k^dkP„,{k)\Fr{k) 



F,^{k)^2 



aHx 



- 1 



F^\k)~{[g + l\D^_),ji{i) 
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+ 1(1 + 1) dx^^D^_ji{x)+2 -^D^_ji{x) 



XX 



dx 
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(B9) 



These expressions are for a sharp source redshift. It is straightforward to generalize them to a distribution dN/dz of 
source redshifts (where dN/dz is assumed normalized), following the discussion in IIV Al In particular, contributions 
to Fj^ (k) of the form 



dxW{x,x)Q{x)ji(x), 



(BIO) 



where Q{x) is any derivative operator, are generalized to 

''''dxW{x,x)Q{x)ji{x)^ 



dxW(x)Q(a;)ji(x) 



pQO ri AT 

mx)= / d~z—W{x = xiS),x)- 



(Bll) 
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Further, contributions evaluated at the source are generahzed to 



Aix)Q{x)ji{x) ^ / dxW{x)Q{x)ji{x) 
Jo 

>V(X)= -^Hiz) 



(B12) 
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